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Abstract 

We discuss many-body states and the algebra of creation and annihilation operators 
for particles obeying exclusion statistics. 
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Recently Haldane introduced a variant of fractional statistics for which the notion 
of exclusion, a generalization of the Pauli principle, rather than exchange or braiding 
properties, is the prime characteristic M. This exclusion statistics appears naturally in 
the fractional quantum Hall effect, spin-| antiferromagnetic chains and the Calogero- 
Sutherland model [1_5 l 

Exclusion principle implies that the number of available one-particle states should 
change with increasing occupation of a state and hence exclusion statistics can be char- 
acterized by the change Adi in available states as the occupation number is changed by 
ANi, i.e., by 

Adi = -J2dij (1) 
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gij are the parameters characterizing the statistics, di is the number of one-particle states 
available to the iVj-th particle with quantum numbers i, holding the labels of the (iVj — 1) 
particles fixed. Particles obeying exclusion statistics were called #-ons in ref.[6], a name 
we shall also use. 

The thermodynamic distributions for particles obeying exclusion statistics can be 
set up generalizing the long familiar combinatorial calculation of entropy for bosons and 
fermions. A number of thermodynamic properties have been studied in this framework! 4-6 ]. 
A question which naturally arises in this context is whether one can go beyond the ther- 
modynamic formulation, attempting a microscopic description by explicitly constructing 
the many-body Hilbert space and further introducing creation and annihilation operators 
and appropriate commutation rules. 

In this paper we discuss an algebra of creation and annihilation operators which 
realizes some of the general features following from Eq.(l). The operator algebra is dif- 
ferent from algebras discussed in connection with the Calogero-Sutherland system^ 7 !; fur- 
ther it seems to be different from proposals for anyon commutation ruW 8 ] as well as 
(/-deformations of boson and fermion algebras, although there are many similarities^ 9 ' 10 ]. 

We start by considering g-ons of a single energy u, i.e., only one species. Let \n) 
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denote an n-body state of g-ons. (We consider them to be in a cavity of volume V = L 3 
and fixed momentum of the form ^ x integer), n can be taken as the eigenvalue of some 
hermitian operator N. Further, we can assume (n\m) = S njm . 
We now introduce an annihilation operator a by 

a\n) = f n \n - 1) (2) 

By orthonormality of states, (n — l\a\n) = / n , which leads, by the definition of the adjoint, 
to 

a t |n) = /: +1 |n+l) (3) 

For a sensible annihilation operator we need a|0) = 0, or f = and the one-particle state 
can be defined as a'\0), so that /j* = 1. 

The number operator N, by definition, must obey the commutation rules 

[N, a] = -a, [N, a f ] = a f (4) 

from which it follows that [N, a^a] = [N, aa^] = 0. The number operator iV can thus be 
written, in general, as a function of a^a or aa) . 

From Eq.(l) we see that for g = 1/m, m being an integer, we lose a state every time 
the occupation number increases by m. We thus consider that in a microscopic description, 
Eq.(l) implies that the maximal occupancy of the state is m. In terms of operators this is 
equivalent to the condition (a^) m+1 = (a) m+1 = 0. A basis for the Hilbert space consists 
of the states |0), 1 1 ),..., |m). By definition the number operator iV satisfies the polynomial 
relation iV(iV — 1)(N — 2)...(N — m) = 0. Further it can be written as an m-th order 
polynomial in a)a (or aa^), 

m 

N = J2M^a) k (5) 
fc=i 

Applying this on |0), |1), |2), etc., and using Eqs.(2) and (3) we get the equations 

m 

J2\fn\ 2k ^k = n, n = l,2,...,m (6) 

fc=i 
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We can now solve these equations for X k and hence obtain N as a function of \fi\ 2 , I = 
2, to. Similarly we can write down a "commutation rule" expressing the relation between 
aa) and a^a, 

m 

aa) = 1 + ^ c r fc(a t a) fc (7) 

where oVs satisfy the equations 



k=i 



^ l/n| 2fc CTfc = |/n+l| 2 - 1 



(8) 



A;=l 



The general solution to Eqs.(6) and (8) is given by 



m 



/ 1 

kldetM 



Qh 



(dxi)- 



-detM 



(xi+i - 1) 1 
fc! detM 



Qh 



(dxi) 



-detM 



x t =0 



where 



detM = —(xiX2—x m ) Y\_( x i ~ x j) 

i<j 



(9) 



(10) 



and Mik = {xi) k . We substitute xi = \fi\ 2 after evaluating the right hand side of Eq.(9). 
For example, for to = 2 ( "semions" ) we find 



aa' 



\h 
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-1/2 
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+ 1 


1/2 


4 -1/2 
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-1/2 


2 


+ 1 


1/2 
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(a f a) 2 

1 - (1 - |/ 2 | 2 ) ata - 1 + {lf f~ 1)1/2,2 (at) 2 a 2 



1/2 



(11) 



iV = 



1/; 


d 4 -2 


l/ 2 l z 


l -|/ 2 | 2 



t , 2-|/ 2 
a'a + 



t ,2-|/ 2 | 
a'a + 



l/ 2 | 4 -|/ 2 | 



t^2 



(a T a) 



1/2 



(a T )^a 



(12) 



We further have (a f ) 3 = 0, a 3 = 0. (The case to = 2 is similar to the case considered in 
ref.[9].) 

Eqs.(5-10) give an (to-1) -parameter set of operator algebras, the parameters being 
x l = \fl\ 2 i I = 2,..., to. Specializations of this class of algebras can be obtained by 
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imposing relations among \fi\ 2 . Such relations can be written in the form of a recursion 
rule 

\f n+1 \ 2 = K(n,g)\f n \ 2 , n=l,... (13) 

the choice of the function K now parametrizing the possible algebras. Once a K(n,g) 
is chosen, Eq.(13) can give, upto a phase, a recursive determination of / n 's and hence a 
representation for a and a*. Some restrictions on K(n,g) follow from general arguments. 
In order to obtain f m+ i = from Eq.(13), consistent with the maximum-m occupancy for 
a state, K(n, g) must vanish for n = m = 1/g. Thus we may expect a factor like (1 — gn) 
to some positive power. Further we have the limiting cases K(n,g) = for g = 1 and 
K(n,g) = for (7 = 0, corresponding to fermions and bosons respectively. |/ n +i| 2 is 
essentially the probability of introducing an extra g-on into a state with occupation number 
n. For m very large and n/m small, the system is very nearly bosonic and the exclusion 
effect is insignificant; we expect |/ n+ i| 2 ~ (1 — gn) + n, corresponding to a spontaneous 
emission term proportional to the available states (1 — gn) and a stimulated emission term 
~ n. Examples of possible choices for K(n,g) are 

g( „, g) = (i-g")[i + (i-gW (14a) 
g( „, g) = (i- g ")'["+(i- g H 1 -' (146) 

There are of course many other choices possible. We expect \f n +i\ 2 /\fn\ 2 to be proportional 
to the number of ways of introducing an extra g-on into the n-filled state, which is related 
to the combinatorial entropy per particle, say S(n) : as ~ e S(n ). The combinatorial rule 
introduced in ref. [4] gives the entropy as 

'l + (l-g)n 



S = log 
S(n) = log 



1-gn 

{l-gny[l + {l-g)nY-v 



+ nS(n) 

(15) 



n 

Eq.(15) suggests the choice (14b). However we would like to emphasize that there is no 
compelling argument for such a choice. In particular the thermodynamic description of 

4 



ref. [4] cannot be derived from a microscopic theory in the standard way, given that there 
is no positive density matrix associated with it^. We cannot therefore expect to recover 
Wu's thermal distribution from Eq.(14b). (Note that the maximum occupancy being m is 
sufficient to ensure that the thermodynamic distribution for the number density n — > m 
as the temperature goes to zero.) Nevertheless we shall analyze choice (14b), and its 
appropriate generalization in the case of many species, in some detail below because it is 
consistent with expected limiting behaviours and interpolates smoothly between bosons 
and fermions. Other choices will give qualitatively similar features. 

A further restriction on K(n,g) is derived from the fact that |/| 2 , being probabilities, 
cannot be negative. A closer inspection of Eq.(14b) reveals that there are inconsistencies 
for some values of the parameter g. The inconsistency appears because, for some values of 
g, one can eventually encounter negative \f\ 2 . Our approach from now on will be that we 
shall use Eq.(14b) only for g = l/m. 

Using the recursion rule (14b) we find 

/i = l 

|/n| 2 =n^ ^ — > 2<n<m (16a) 

k=i 

= 0, n > m + 1 (166) 

As we mentioned earlier, this provides, upto a phase, a matrix realization of a and at; 

*mn $m,n — 1 fni ^rnn $m,n+l fn+1' 

We now turn to the case of many species corresponding to different values of energy or 
other quantum numbers. For two species of g-ons we introduce creation and annihilation 
operators a\,ai, i = 1,2 and fn},n 2 with 

a 1 \n 1 ,n 2 ) = fi\] n2 \ni - l,n 2 ), 

a 2 \n l ,n 2 ) = f n % 2 \n 1 ,n 2 - 1), (17) 
where f^l 2 = f^o = an d /i*d = fo^i = 1- Eq.(3) can be generalized to the case of 



two-species as 



\f { n^i,nf = K^\n u n 2 ,g)\f^ n2 \\ m >0 

\tiL+i\ 2 = ^ (2) K,n 2 ,^)|/S n2 | 2 , n 2 > (18) 



A suitable generalization of Eq.(14b) is 



K^(n u n 2 ,g) = (l + w i )l[(- 



\9ki 



k - +Wk 



WiUi + gijUj = 1 (19) 

The many-body states and the representation of and a\ can be constructed by solving 
these recursion relations. We shall first address the case where gij is diagonal, i.e., no 
mutual statistics, and for the parameters gu = 1/mwe shall obtain the full operator alge- 
bra. We shall later discuss the case with mutual statistics and the somewhat complicated 
structure of the corresponding many-body Hilbert space. 

In the absence of mutual statistics, the recursion rules (18) reduce to 

(1) . |2 = ai h ,„ u \l-gilf(l) |2 



n 2 I 



l/Sn 2 + l| 2 =<(l + ^2) 1 --|/S n2 | 2 (20) 

where gi = gu and riiWi = 1 — girii. For consistency reasons outlined earlier in the 
case of one species of (7-ons, we shall consider special values of g^ in particular the case 
9i =92 = 1/m. We again find that the maximal occupancy for each species is m; i.e., 
(ai) m+1 = (aj) m+1 = 0. Further, the recursion rules (20) are sufficient to show that the 
Hilbert space is a product of the Hilbert spaces for each species. These properties clearly 
generalize to arbitrary number of species. 

The recursion rules (20) do not completely determine the functions fn},n 2 - Specifically, 
l/i ^J and l/^il for ni, n 2 7^ 0, as well as the phases of all the /'s, are undetermined. 
Thus Eq.(20) does not suffice to get all the matrix elements of ai and a\. This is, of 



course, no different from the situation with bosons or fermions; one needs commutation 
rules relating operators for different species in order to obtain a complete operator algebra 
and determine all the matrix elements. 

The rationale for the choice of off-diagonal commutation rules has to come from the 
exclusion principle in the following way. Although not specified explicitly, we have been 
dealing all along with free g-ons. Thus the quantum number labelling the species is the 
spatial momentum p with a corresponding energy u(p). The property of exclusion reads, 
for <7j = 1/m, (a p ) m+1 = (aj } ) m+1 = 0. Eventhough the p-diagonal representation is 
appropriate for free particles, we can equally well use the coordinate representation. Since 
the coordinates x also provide a complete set of quantum numbers for the particles, we 
expect the exclusion property (a(x)) rn+1 = 0, (a^(x)) rn+1 = to hold, where a(x) = 
e ipx a p and a'(x) is the adjoint of a(x). Of course, one could also use a(q) = J2 q u q (p)a p 
in (/-diagonal representation, q's being any complete set of quantum numbers and u q (p) 
being the appropriate functions. Thus in general we expect linear combinations of a's to 
obey the exclusion principle, i.e., 

(^a^r +1 = E<aJr +1 = (21) 

i i 

where are arbitrary. 

From the definition of the number operator Ni\n\, ri2, ■■■) = ni\ni, ^2, •••)> we find that 

[N i ,a j ] = [N i ,a]} = 0, i^j (22) 

Eqs. (21,22) will be our main guide for obtaining the full operator algebra for the case of 
many species. However they are still not sufficient to determine the algebra completely, so 
we shall further assume that 

didj = Rijdjdi, i ^ j (23) 

where Rij is a c-number. Consistency of Eq.(23) requires RijRji = 1 . Further, d\d - = 
Rji a j a t- Since iVj is a function of d\d i: we seek to satisfy Eq.(22) with the slightly stronger 
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requirement [aja^, aj] = [ajaj, aj] = 0, z 7^ j. The general solution to this equation, using 
Eq.(23), is 

a l a j = R i 3 la A 

cua) = Rijd^ai, i ^ j (24) 

Consistency of these equations gives R*jRij = 1- Thus Rij is a phase, R^ = e %dii . We 
must now impose Eq.(21). For two species, we have 

m 

(a 1 a 1 +a 2 a 2 ) rn+l = {a 1 a 1 ) m+1 + (a 2 a 2 ) m+1 +Ys C ( k ' m ' ^ 2 )((«iai) m+1 " fc («2a 2 ) fc (25) 

k=i 

where 

In particular C(l,m,i2) = Since (ai) m+1 = (a 2 ) m+1 = 0, Eq.(21) can be 

satisfied if all the C(k,m, R* 2 ) vanish. This determines Ri 2 . The generic solution for 
arbitrary m is of the form R\ 2 = exp(±2-7ri/(m + 1)). There can be more solutions for 
special choices of m. All solutions are a subset of the (m+l)-th roots of unity. For example, 
for m = 2, R 12 = e ±27 "/ 3 ; m = 3, R 12 = e ±2 ™/ 4 ; m = 4, i? 12 = e ±2iri / 5 , e ±4 ™/ 5 . Since 
we want to treat all m's on an equal footing, we are going to use the generic solution 
-R12 = exp(±2-7rz/(m + 1)). The commutation rule becomes aia 2 = e 27Tl ^ m+1 ^a 2 ai or 

Notice that the second choice is obtained from the first by a 
relabelling a\ <-> a 2 . The basic solution will be R\ 2 = e 27r V( m + 1 ) ? U p to such relabellings. 

When we have more than two species, the exclusion condition (21) can be satisfied as 
follows. Write 

aia 3 = Ri 3 a 3 ai, a 2 a 3 = R 23 a 3 a 2 (27) 

For the choice R 23 = -R13, 

Aa 3 = R± 3 a 3 A 



where A = a\a\-\-a 2 a 2 . The binary expansion (25) suffices to simplify (A + a 3 a 3 ) m+1 = 0. 
We get the same solution as before, viz., i? 13 = i? 12 ; the choice Ri 2 = R13 = R23 obviously 
generalizes inductively to arbitrary number of species. 

Recapitulating, the commutation rules we have so far are, for Oj = 1/m, 

ai a\ = 1 + ^ (y k {a\a,) k (28a) 
k=i 
2ni 

CLidj = exp( — — ^-J ajdi 1 < j (286) 

a\a] = exp(-^-) a)a\ i < j (28c) 
+ . 2ni . + . . . 

a\aj = ex p(-^qrr) a J a l 1 < j ( 28d ) 

We have arrived at Eqs.(28) by using the simplifying ansatz (23), with Rij being a c- 
number. In order to write these commutation relations for different species we have to 
introduce an ordering, for example i < j. The indices etc., being momentum labels, 
ordering is naturally possible only in one spatial dimension. Since many of the physical 
situations where exclusion statistics might be relevant are effectively one-dimensional^ 1-5 ], 
this may not be a drastic limitation. It is possible that more general structures can be 
constructed taking R^ to be an operator. Notice also that once we have chosen an ordering 
of momenta we cannot relabel a± <-> a 2 , etc. The commutation rules with R^ <-> R*- are 
thus distinct. 

(i) 

We can now write down /Ai,n 2 thus providing an explicit matrix realization for a\ and 
a\. Consider two species, i = 1,2. From Eqs. (17,28) we derive two sets of equations that 
should be satisfied by /'s, 

fi%A\] n2 -i = exp(27u/(m + l))/ff n2 /£ } _ 1)n2 (29a) 
f&f&r* = exp(-27rz/(m + l))fV n2 _J* n %, n2 (296) 
Eqs. (20,29) completely determine the moduli |/|; in particular, 

1/SnJ 2 = Ifi'L-il 2 = - = l/Sol 2 = l/nj 2 , < m.na < m 
1/SnJ 2 = l/i 2) -i,nj 2 = - = l/o ( ll 2 = l/n 2 | 2 , < n u n 2 < m (30) 
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where \ f ni \ 2 are defined in Eq.(16) for a single species. 

Eqs. (29) also determine some of the relative phases of the /'s, but they still allow 
some freedom of choice for the phases; this is related to the ambiguity in how many-particle 
states are defined. For example, (a|) ni (aj) n2 10) and (4) na (a{) ni |0) have the same number 
of particles but different phases. We will make the choice 

|n 1 ,n 2 ,...> = C(aJ) ni (4) na -|0> (31) 

where C is a normalization factor chosen to be real. This is ordered in terms of ascending 
momenta; the a^'s with smaller labels appear to the left. With this phase convention, we 
have, for two species and m = 2, 

f (i) _ f« _ f (i) _ f (2) _ i 

J 1,0 ~ J 1,1 ~ J 1,2 ~ J 0,1 ~ L 

f(l) _ fU) _ f(l) _ f(2) _ 
/2,0 ~~ /2,1 — /2,2 — J0,2 — 

f (2) _ ^27ri/3 f (2) _ 2iri/3 

J 1,1 ~ e ) Jl,2 — e 



1/2 

(32) 

All other /'s are zero. We thus have a complete specification of a i7 a\. This can be 
easily extended to many species. In fact, it is evident that the commutation rules (28) 
suffice to evaluate the action of a^s and aj's on any state of the form (31). The exchange 
properties of the many-body wavefunctions are evident from the representation (31) and 
the commutation rules (28). 

In the case of mutual statistics, viz., g^i not diagonal in Eq.(19), the structure of 
the many-body Hilbert space is somewhat more complicated. The number of many-body 
states can be greater or smaller, depending on the signs and magnitudes of g^i, than what 
is expected from the tensor product of states for each species. This change in the expected 
number of states has to do with Eqs. (18, 19). From Eq.(19), Wi can become negative 
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r( 2 ) _ g47ri/3 



J2A — 



f(2) _ 
J2,2 



as rij increases. Since \ f^\ 2 is positive, a negative Wi can give an inconsistency, as we 
have already discussed earlier, unless some of the /'s are zero. Notice that Eq.(18) leaves 
/i,n 2 ' fnuii f° r n i' n 2 7^ undetermined. One can set some of these to zero to obtain a 
consistent solution to Eqs.(18,19). The possible states are then no longer what is expected 
from the tensor product. 

We shall illustrate this by considering two characteristic examples with two species, 
where the number of the states increases and decreases respectively, compared to the tensor 
product of the single species states. 

A) 011 = #22 = ^, 912 = 921 = 

In this case the recursion rules reduce to 

j_ 

fii 1 J 2 = (l + ^) 1 -»(^) B (l + ^)»l«J J (33a) 

l^n 2+1 | 2 = (l + - 2 ) 1 -"(^) m (l + -l)"l^n 2 | 2 (336) 

where wini = 1 — ^ ni ~ n2 - > , w 2 n 2 = 1 + ( ' ni ~ n2 ' > . We see that the ratio w\jw 2 becomes 
negative if \n\ — n 2 \ > m. A consistent solution requires that fim+i = and f^+i 1 = 
0, i > 1. This choice allows only nonnegative wi/w 2 . (This also avoids the singular points 
«)2 = in Eq.(33a) and w\ = in Eq.(33b)). The allowed states are of the form \ni,n 2 ), 
where |ni — n 2 \ < m. There are (2m + l)(m + 1) such states, which are m{m + 1) more 
than in the tensor product of (m + 1) states for each species. 

B) #11 = #22 = ^, 912 = 921 = ^- 

In this case the recursion rules become 

\C } + i,nf = (l + ^i) 1 "- (wiw 2 )^ (l + w 2 )-^\fi\] n f (34a) 
\fi 2 L+i\ 2 = (l + w,) 1 -^ (w lW2 )^ (l + w 1 )-^\fi%f (346) 

where wini = 102^2 = 1 — ^ Ul ^ 12 ^ ■ Here we must require /i^+j = and f^+i 1 = 0, % > 0, 
so as to avoid recursive sequences which lead to the singular points w 2 = — 1 for Eq.(34a) 
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and wi = — 1 for Eq.(34b). The allowed states now have (ni + n 2 ) < m. There are 
\{m + 2)(m + 1) such states, which are \m{m + 1) less than in the tensor product of 
(m + 1) states for each species. 

It is clear from the above discussion that the operator algebra for the case of mutual 
statistics must be fairly involved so as to reflect the complicated structure of the many- 
body states. This is currently under investigation. 

Acknowledgements: The commutation rules (5,7) can be mapped into the usual 
bosonic algebra^ 11 ]. We thank A. P. Polychronakos for discussions and for pointing this 
out. 

This work was supported in part by the DOE grant DE-FG02-90ER40542 and the 
NSF grant PHY-9322591. 

Note added : The combinatorial argument after Eq.(14), with equal a priori probability 
and maximal occupancy m, gives yet another choice for K{n,g = 1/m) as 

1 m+ 1 

K - 1 if ™+! - 1 

It is unclear how to generalize this equation to the case of mutual statistics. The case 
m = 2 for the above choice of K is similar to the case considered in ref.[12]. 
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